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One of the most frustrating issues in early universe cosmology centers on how to reconcile the vast 
choice of universes in string theory and in its most plausible high energy sibling, eternal inflation, 
that jointly generate the string landscape with the fine-tuned and hence relatively small number of 
universes that have undergone a large expansion and can accommodate observers and, in particular, 
galaxies. We show that such observations are highly favored for any system whereby physical 
parameters are distributed at a high energy scale, due to the conservation of the Liouville measure 
and the gauge nature of volume, asymptotically approaching a period of large isotropic expansion 
characterised by w = — 1. Our interpretation predicts that all observational probes for deviations 
from u) = — 1 in the foreseeable future are doomed to failure. The purpose of this paper is not 
to introduce a new measure for the multiverse, but rather to show how what is perhaps the most 
natural and well known measure, volume weighting, arises as a consequence of the conservation of 
the Liouville measure on phase space during the classical slow-roll expansion. 

PACS numbers: 04.60.Pp, 98.80.Cq, 98.80.Qc 


Cosmology in its standard form, the so-called AC DM 
model, provides a remarkably successful fit to all cur¬ 
rent data. The latest developments, from the large-scale 
structure of the galaxy distribution, dark matter sur¬ 
veys via gravitational lensing, dark energy via supernovae 
and cosmic microwave background fluctuations via or¬ 
bital and suborbital telescopes, have ushered in the age 
of precision cosmology. Inflationary cosmology has been 
especially adept in pushing the frontier of our knowledge 
of the initial conditions back to the epoch of quantum 
gravity, the Planck epoch, where one enters the intellec¬ 
tually appealing domain of string theory. However this 
has occurred at a severe cost, namely that of introduc¬ 
ing the landscape of some 10^™ Calabi-Yau manifolds, 
any one of which is equally likely to contain our present 
universe in the absence of any selection principle. This 
is a potential disaster, since galaxies could not form in 
almost all of these future universes, and has led to the in¬ 
troduction of anthropic reasoning as a means of selecting 
a universe with cosmological constant comparable to the 
observed value [l[, observed to be exceedingly small in 
terms of fundamental units from a set where the typical 
value is larger by 120 orders of magnitude. 

For many physicists, this is not a satisfactory reso¬ 
lution. More significantly, in the absence of a measure 
that enables us to assess the probability of any given 
universe in the multiverse, the problem is not even well 
posed. There have been numerous attempts to define 
the measure of the universe, but to date, none have been 
completely successful. Moreover, many of these measure 
prescriptions rely on implementation via anthropic rea¬ 


soning. This is equally true in string theory 13,131 , eternal 
inflation and loop quantum gravity f7Hloj|^ 

Here we propose a novel solution that accounts for the 
size and even acceleration of the universe. Our idea cen¬ 
ters on volume-weighting, and has been discussed previ¬ 
ously iEH- However we show here that it arises natu¬ 
rally as a result of a dynamical systems approach to cos¬ 
mology. Our goal is to show that the volume weighting is 
a direct consequence of the conservation of the symplec- 
tic structure under Hamiltonian flow during the classical 
expansion of the universe. Volume weighting has been 
examined extensively in the context of eternal inflation, 
sometimes referred to as a ‘Stationary Measure’. In 121 
it was shown that, in the case of simple models at least, 
such a measure produced gauge invariant results - a key 
issue plaguing many other proposals such as time or scale 
factor cut-offs. In particular, given a set of possible de- 
Sitter vacua, the probability of an observer finding them¬ 
selves in such a vacuum is proportional to the amount of 
volume growth during the slow-roll expansion of the vac¬ 
uum. In [1^ this proposal was shown to be free of both 
the youngness paradox (that is there is not a high like¬ 
lihood that our bubble was very recently created) and 
of Boltzmann brains. The analysis we present here is 
complementary to such approaches; from the perspective 
of eternal inflation, we will examine the behaviour of a 
single bubble and the evolution of a measure during the 
slow-roll inflationary phase. 


The physical observables defining the universe are de¬ 
termined according to some unknown probability distri¬ 
bution at an initial energy density. This is motivated 
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by the fact that general relativity coupled to the stan¬ 
dard model is an excellent fit to observational evidence 
at low energies in the classical regime. We consider the 
probability of an observer at low energy density making 
physical measurements which are compatible with a given 
evolution of the universe. The strategy is the following: 
we express general relativity (GR) coupled to matter in a 
Hamiltonian setting and separate phase-space variables 
into those which form physical observables and those 
which are gauge. We take the probability distribution 
on physical observables, expressed in terms of non-gauge 
variables, and ‘thicken’ this to a form which is propor¬ 
tional to the Liouville measure. We do not base our mea¬ 
sure on the Liouville measure itself, but rather observe 
that any measure formed from physical observables can 
be thickened to one proportional to the Liouville mea¬ 
sure by taking exterior products with gauge directions to 
cover the whole of phase space. This measure is then con¬ 
served through evolution, and we can thus project back 
down by integrating out the gauge directions to form a 
new probability distribution at late times to describe the 
observations made at low energy. Thus we express the 
induced low density probability distribution in terms of 
the high density distribution, using Liouville’s theorem 
to translate between the two. Critically, the length of 
the interval in gauge direction may not be conserved un¬ 
der evolution, thus on projecting back to physical observ¬ 
ables the interval in this gauge direction must be adapted 
to cover the same set of solutions as chosen at the high 
energy surface. This induces a weighting on the low en¬ 
ergy probability distribution proportional to the relative 
expansion of gauge directions. 

In Q, the Liouville measure was directly applied at a 
low energy density to claim that inflation was an improb¬ 
able event. The discrepancies between this and the cor¬ 
responding result for high energy was explored in [ 3 , 0 | ; 
and the resulting interpretation as an implicit probabil¬ 
ity weighting given in [9|. Our approach differs in that 
we allow for any choice of measure at high energy and 
show that when translated to observations at low energy 
a re-weighting based on increase in volume between equal 
energy slices arises naturally. This differs significantly 
from volume-weighting directly which is introduced 
by hand. 

We will consider an action consisting of matter satisfy¬ 
ing a Lagrangian Cmatter minimally coupled to gravity on 
a topology that splits into a spatial slice E and a time-like 
direction. The matter Lagrangian will be kept general - 
we shall only require that it obeys the weak energy con¬ 
ditions. A more inflationary focussed exploration of such 
dynamical systems is given in 1^ . For brevity of exposi¬ 
tion, we will take gravity to be described by GR, however 
a wide class of gravitational theories exhibits qualitita- 
tively similar properties , and we shall assume E to 
be flat and compact, although these assumptions can be 
relaxed in appropriate limits. 


S = jdt 

matter ) ( 1 ) 

We shall assume that our manifold admits a constant 
mean curvature foliation [l6| . This is not highly re¬ 
strictive, and largely consists of assuming weak energy 
conditions on our matter Lagrangian. Thus the mean 
curvature (Hubble parameter) H will be monotonically 
non-increasing between slices and we can reparametrize 
our time interval as an interval in H. Since we are con¬ 
cerned with properties observed at a given energy density, 
or equivalently, Hubble parameter, we identify a dynam¬ 
ical trajectory as a single entity which is determined by 
the values of physical parameters at a given iJ, from the 
set of all trajectories which obey the Hamiltonian con¬ 
straint. 

This choice of time foliation by Hubble parameter is 
presented for reasons of brevity, and is not strictly nec¬ 
essary. From a physical perspective, any monotonic pa¬ 
rameter can be chosen to form a good clock at least in the 
classical regime. During the highly quantum mechanical 
phase of eternal inflation, the Hubble parameter fails to 
satisfy this condition, but a different choice of time folia¬ 
tion can be employed - the standard route being the use 
of a massless scalar field minimally coupled to gravity. 
Following such a procedure produces an analogous result 
to that given here. 

The action carries naturally a two-form on E, the sym- 
plectic structure w, which is conserved between Gauchy 
slices. We further may decompose our fields into eigen¬ 
functions Ei of the Laplace operator compatible with E 
(Fourier modes) with amplitudes Ai. To make contact 
with observations, we will restrict ourselves to a hnite 
number n of such modes across all fields. Thus our sym- 
plectic structure can be expressed in terms of functions 
fj of the amplitudes: 


u; = Y^fj {A)dTrj A dAj (2) 

3 

where tt^ is the conjugate momentum to Aj. We form 
the Liouville measure by raising the symplectic structure 
to a high enough power to cover the space of physical 
trajectories on a given Hubble slice: 


D = w" (3) 

Of particular importance for our analysis is the back¬ 
ground volume mode u, whose contribution to the sym¬ 
plectic structure is dvAdH. This mode itself is not an ob¬ 
servable, but rather a gauge choice. In order to measure 
any physical length scale within cosmology, a reference 
scale must be given. A ‘rods and clocks’ construction of 
cosmology must concede that the rods and clocks used 
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are present in the evolving universe. Therefore there will 
exist a gauge redundancy under v —>■ \v. One may choose 
a length scale in terms of the Fourier modes, but all that 
can be observed is the ratio of the wavelength of the mode 
to the total volume v, and thus an indistinguishable de¬ 
scription of the same physical system will exist with a 
rescaled volume. It should be unsurprising to learn that 
the same physical system can be described in either met¬ 
ric or imperial units - although a physical length such 
as the volume of a closed universe will remain a physical 
observable, the choice of rod by which it is measured, be 
it the metre or foot, is a redundancy. This is reflected in 
the Hamiltonain constraint which can be expressed 

0 = 'H = vF[H,Aj,Aj] (4) 

and all physical observables at a given Hubble slice are 
composed of Aj , Aj. Since a factor of v can be pulled out 

of all observables, we find tt^ = vfj[A,A] for some func¬ 
tions fj. We require that the Hamiltonian constraint is 
satisfied, and so not all the Aj and iTj are independent. 
We can therefore eliminate one of these using the con¬ 
straint, and thus (up to an arbitrary choice of sign) the 
Liouville measure becomes 

H = Aodv A dAo /y diTi A dAi (5) 

i^o 

and we can use the constraint to express Ao in terms of 
the remaining modes and the Hubble parameter. Rewrit¬ 
ing this measure in terms of the 2n — 1-form $ on the 
remaining phase space directions, we can express the Li¬ 
ouville measure as 

= M[H, A, A]d^ A dv (6) 

wherein M is the expression of Ag in terms of the re¬ 
maining physical observables. Let us now construct the 
probabiltiy of observing an event £. We suppose that the 
physical parameters in the early universe are distributed 
on an initial Hubble slice Hi according to some unknown 
probability distribution p, thus the probability of observ¬ 
ing £ is 

P{£)= [ p[H,Aj]d<^> (7) 

JR 

in which R is the region of physical parameters for which 
the event £ occurs. We are free to ‘thicken’ this region 
by an interval / in volume, which we shall take to be of 
unit length, to obtain an integral that is proportional to 
the Liouville measure (dropping the parameters from p 
and M): 

P{£) = j pd<^>Adv= f (8) 

H R'X, I J R X I 

An observer wanting to recreate this probability by mak¬ 
ing observations at a later time, and thus lower mean 


curvature H' can do so defining R’ to be the region to 
which trajectories in R have evolved: 


P{£)= [ p'[H',A,A]d^ (9) 

JR' 

To obtain p' from p, note that Liouville measure is con¬ 
served, therefore we can again integrate this over a range 
I' of the volume direction to eliminate the gauge direc¬ 
tion, and recover a measure on the physical parameters. 
However, although v is gauge in terms of observables, the 
relative expansion in the gauge direction is dependent on 
the physical parameters: To cover the same set as trajec¬ 
tories as was in i? X / at the initial surface, we must allow 
I' to vary. The trajectories lying in the interval [ni,r! 2 ] 
on Hi now occupy the range [Afi,Ar! 2 ] with A the ratio 
of volumes between slices: 


A = 



HdH 

H 



dH 

H{l + w[H]y 


( 10 ) 


in which we have expressed the ratio between the homo¬ 
geneous components of pressure and energy density as 
w = P /for the effective Friedmann equation describ¬ 
ing the evolution of the homogeneous modes 2^. In eq 
(fTOl) the first line is generic across all gravitational the¬ 
ories as it represents a geometric identity, whereas the 
second represents the application of GR. The precise for¬ 
mulation of w[i7] is dependent upon the matter model 
and initial conditions, but this formulation holds for all 
weakly coupled matter models obeying the weak energy 
condition. Thus to recover P{£), we must integrate the 
Liouville measure over a range in volume. Since the 2n- 
form separates into dv A d^, the volume integral can be 
performed hrst. Thus to obtain the probability of find¬ 
ing physical observables in state X at Hi, which evolve 
to X' and H 2 according to a probability distribution at 
high energy density, we must multiply the low energy 
distribution by the volume expansion: 


pAx) = x{xAp{x') (11) 


and so the probability must be weighted by the expansion 
in volume observed. It is easy to see from eq. m that 
those solutions which have w ~ — 1 for a large period of 
their evolution must be highest weighted, and that this 
weighting can be very large. In toy inflationary mod¬ 
els consisting of a scalar field with a quadratic potential 
(m = 1.21 X 10~^Mpi), taking the initial Hubble slice to 
be at the Planck scale this ratio reaches exp[10^^]. Thus 
the weighting can be completely overwhelming relative 
to any other measure considerations. If one pushes the 
initial distribution all the way back to the singularity in 
GR, to recover the weighting becomes a delta function on 
those trajectories which have maximal expansion. This is 
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of course a purely mathematical exercise, as dimensional 
analysis informs us that quantum gravity effects must 
occur around the Planck regime, and thus our classical 
model is invalid here. 

Anisotropies contribute a term that falls off as 
(equivalently for scale factor). This can be seen 
directly from the Bianchi I metric which describes the 
background modes of an anisotropic, homogeneous, flat 
cosmology in terms of the anisotropy parameters j3 and 
cr. The spatial part of the metric has line element: 

dsl = ( 12 ) 

When inserted into our action O/? and <7 have the same 
form as massless scalar fields. Note that in this action, 
the values of a and /3 are gauge, as a rescaling of coordi¬ 
nates can change them arbitrarily, however their kinetic 
terms are observable and contribute to dynamics. Again 
we must integrate the Liouville measure over a range in 
these gauge directions to project back down to a measure 
on physical observables. In this case, however, the inter¬ 
val is preserved under evolution, and so this is a trivial 
move. Those solutions which have high values of /3 and & 
undergo far less expansion than more isotropic solutions, 
and are therefore disfavored after the effect of expansion 
in the volume gauge direction is considered. This pro¬ 
vides a counter to the claim that isotropic universes are 
a set of measure zero - as one takes the infinite limit 
of initial energy slice the distribution at high energy be¬ 
comes a delta function on isotropy. 

A similar argument holds for inhomogeneities: Inho¬ 
mogeneities distributed (as Fourier modes) at an initially 
high energy slice would lead to less expansion between 
surfaces of constant mean curvature. Therefore on ex¬ 
amination of the distribution of anisotropies at low den¬ 
sity, smaller amplitudes of inhomogeneous modes would 
be more strongly weighted to recover the high density 
distribution. This can be seen by considering the space 
of perturbative modes about a homogeneous (FRW or 
Bianchi I) background, but such an analysis ignores back- 
reaction. The result holds for a constant mean curvature 
slicing in generality, and is therefore far stronger. This 
can be further generalized by substituting any monotonic 
observable (e.g. massless scalar field) to play the role of 
clock in situations where the mean curvature is not suit¬ 
able. 

Volume weighted measures are often criticized for en¬ 
countering the “Q-catastrophe”, in which it is noted that 
maximal expansion disfavours the amplitude of density 
fluctuations, Q, being the observed 10“®[l^. This value 


assume an approximately flat prior on the mass (or other 
coupling constant) of the scalar field which is assumed to 
be the inflaton, and analyse the resulting the resulting 
distribution of Q, leading to a strong preference for large 
values of Q. However, it has been noted that a flat prior 
is not justified for this parameter as zero is a special 
point of the distribution, and in , a more refined dis¬ 
tribution is explored. Here it is noted that the problem is 
inverted - such an approach leads to lower values of Q be¬ 
ing preferred, with I0“® being orders of magnitude more 
likely than the observed 10“^. Such issues are avoided 
in the multiverse context by the application of a scale 
factor cut-off 2l| , however such a cut-off is incompatible 
with evolution of the Liouville measure - different initial 
conditions will cross any cut-off at a rate dependent on 
their physical parameters. A separate solution is offered 
through decay of the inflaton [23 or through a specific 
choice of landscape. 

In such analyses, it is assumed that inflation is due to 
a scalar field subject to a power-law potential, and that 
this field is also responsible for the perturbation spectrum 
that is observed. There exist a whole plethora of possi¬ 
ble inflationary mechanisms, and the set in which this 
is due to a scalar field subject to power law potential 
are but a small fraction of these. Furthermore, highly 


successful models such as the curvaton [23| and multi- 
field inflationary models 2^ are claimed not to suffer 
from this issue at all: In such models the amplitude of 
fluctuations is decoupled from the behaviour of the infla¬ 
ton, removing the Q-catastrophe. However such models 
require that the mass of the curvaton, for example, is 
forced far below that of the inflaton, and thus the decou¬ 
pling may not be complete [ 2 ^. We note at this point 
that the Planck data indicate that such single field mod¬ 
els are observationally disfavored 23 , and therefore it is 
natural to consider modified scenarios such as multifield 
models. The question of the Q-catastrophe, and its pos¬ 
sible resolution, rests on probability distributions on the 
landscape, which we do not address. The analysis that 
we present here makes no claim on the nature of the in¬ 
flaton - the re-weighting of measures by volume due to 
the conservation of the Liouville form is independent of 
the matter playing the role of inflaton. This can even 
be extended to modified gravitational actions in which 
the inflationary behaviour is due to higher order terms 
of the Ricci scalar (for example) for which there is no 
Q-catastrophe, as all that is required is the conservation 
of the symplectic form, a standard result of Hamlitonian 
physics. 


is in the region considered to be necessary for the creation 
of galaxies in the era before dark energy domination, and 
is correlated with the flatness of any inflationary poten¬ 
tial, with the formation of galaxies stopping if Q < I0“®, 
and failing due to over density if Q > 10“"^ [l9|. The 
standard analysis which leads to the ‘catastrophe’ is to 


The idea of generality is of course vague [26j and of¬ 
ten subject to the biases of the observer’s preferred set 
of variables 27|: One cannot define a probability distri¬ 
bution at high energy without reference to some more 
fundamental theory. There are several alternative ap¬ 
proaches to defining cosmological probabilites, particu¬ 
larly in the areas of inflation. Most of these, however. 
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start from a model which involves directly the gauge di¬ 
rection in volume, leading to a host of measure issues Q. 
These arise most acutely in the context of eternal infla¬ 
tion, in which a (countably) infinite number of “bubble 
universes” can form, and thus probability definitions un¬ 
avoidably encounter classical problems of counting infi¬ 
nite sets in which any probability that is desired can be 
obtained by re-ordering the set. Attempts to solve this 
often involve ad-hoc arguments, such as cut-offs in time 
™ or those involving maximizing numbers of observers 


However, these arguments all have significant flaws: 
There is no way for an observer in a post-eternal infla¬ 
tionary universe to know how long the universe spent 
in this phase, and indeed if the nucleation rate is high 
enough, for any given time T the probability that we 
spent longer than T in this state approaches 1 as the 
cut-off time tends to infinity. Maximizing observers is 
also heavily dependent on the order of counting - it is 
ill-defined to ask whether it is more probable that we 
belong to one of the infinite number of universes with a 
single observer, or a single universe with an infinite num¬ 
ber? The argument put forth here is that one should be 
agnostic on such issues, as they represent unobservable 
parameters. The result presented is that the manner in 
which initial parameters for a post-quantum inflationary 
phase are distributed is largely irrelevant. What matters 
is not what die is cast, nor how it is cast, but rather 
when: by whatever means the physical parameters are 
distributed at some initial energy density, the probabil¬ 
ity that an observer at a lower energy density observes 
an event is weighted by the relative expansion in vol¬ 
ume of the set of universes for which that event occurs. 
This weighting favours heavily those solutions which un¬ 
derwent a large expansion, i.e. for matter obeying the 
weak energy condition, those which had an extended pe¬ 
riod in which w = —1, and are closest to isotropy and 
homogeneity, as these experience the greatest expansion. 
Note that the purely isotropic, homogeneous, inflation- 
dominated universe is a set of measure zero on the space 
of all possible trajectories since these form a continuum, 
so a typical observer should expect to find conditions in 
a small neighborhood around this point, just as an expo¬ 
nential distribution is peaked at zero - this is the most 
probable state, but any measurement should be expected 
to be in a region close to this. 

This inverts the usual claim of fine-tuning of initial 
parameters: it would in fact require that at the point 
of inception, the probability distribution of physical pa¬ 
rameters would have to be heavily fine-tuned to avoid a 
large, isotropic universe which experienced slow-roll in¬ 
flation. Maximal expansion can of course be generated 
by a pure cosmological constant. However, since we are 
interested in the observations made at low energy den¬ 
sity (such as in the universe that we currently inhabit) 
this precludes the case of pA exceeding this value. Such 


reasoning is, of course, anthropic and follows the argu¬ 
ments of Weinberg [l[ - volume weighting of probabilities 
does not explain why the cosmological constant is within 
any set of bounds, but strongly favours the maximiza¬ 
tion of the cosmological constant subject to this condi¬ 
tion. Thus our analysis does not explain the existence of 
such a bound, but rather its saturation. A low energy 
observer should expect to see a universe dominated by 
conditions very close to w = — 1 with H = 1. 

Our derivation of volume weighting applies during the 
classical, slow-roll expansion of the universe regardless of 
any eternal inflationary phase. This derivation comple¬ 
ments the arguments of Linde et al. 12, [^| who pro¬ 
posed volume weighting during slow roll in the context 
of eternal inflation. We have shown that distributions at 
high energy density are re-weighted by the flow of the Li- 
ouville measure under Hamiltonian evolution. Given an 
distribution of physical parameters from some pre slow- 
roll phase (be it based on eternal inflation, bouncing cos¬ 
mologies or any other considerations) our analysis shows 
that on making low energy observations after slow-roll, 
the conservation of the symplectic structure introduces 
naturally a volume weighting. 
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